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Abstract 


A  rigorous  theory  of  backscatteririg  from  slightly  rough  statistically  nonstationary  sur¬ 
faces  is  developed.  The  development  precedes  through  the  application  of  two  ensemble 
averages;  the  first  to  all  realizations  of  the  surface  having  fixed  nonstationary  features,  and 
the  second  over  all  nonstationary  surface  characteristics.  The  resulting  theory  is  applied  to 
acoustic  seafloor  reverberation,  where  it  is  found  that  nonstationarity  must  be  assumed  to 
obtain  agreement  with  observed  reverberation  fluctuations.  The  analysis  is  further  extended 
to  the  time  domain  where  it  is  demonstrated  that  for  Bragg  scatter,  no  reduction  in  the 
predicted  fluctuations  occurs  for  frequency  diverse  waveforms,  thus  providing  a  method  of 
differentiating  Bragg  scatter  from  other  nonresonant  scattering  processes. 
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1.  Introduction 

The  scintillations  of  signals  propagated  through  turbulence  or  scattered  from  rough  sur¬ 
faces,  such  as  the  twinkling  of  stars  or  the  play  of  sunlight  on  water,  have  long  fascinated 
man.  As  a  consequence,  extensive  research  has  been  devoted  to  first  understanding  and 
explaining  such  phenomena  and  subsequently  to  exploiting  them  to  probe  the  environment. 
Despite  these  efforts,  the  nature  of  the  reverberant  field  backscattered  from  rough  surfaces 
at  low  grazing  angles  is  still  poorly  understood. 

An  analogy  is  often  drawn  between  rough  surface  scatter  and  scattering  from  a  random 
distribution  of  discrete  points  [Goodman  1976].  In  this  case,  the  scattered  field  is  regarded  as 
the  complex  sum  of  the  contributions  from  the  many  elementary  areas  of  the  surface  that  are 
simultaneously  irradiated.  If  the  contributions  from  each  scattering  element  are  independent 
and  their  phase  is  uniformly  distributed  over  its  full  range,  then,  as  the  number  of  such 
areas  becomes  large  the  distribution  of  the  amplitude  of  the  scattered  field  approaches  the 
Rayleigh  distribution.  In  this  article,  the  scintillation  index,  which  is  unity  for  the  Rayleigh 
distribution,  is  used  as  a  measure  of  the  strength  of  the  fluctuations. 

Measured  statistics  of  microwave  sea  surface  clutter  and  of  acoustic  seabed  reverberation 
are  generally  found  to  exhibit  significantly  higher  fluctuations  than  those  predicted  by  the 
Rayleigh  distribution.  A  number  of  approaches  have  been  taken  to  model  or  explain  these 
phenomena.  Jakeman  and  Pusey  [1976]  and  Oliver  [1983]  propose  models  that  basically 
stem  from  a  consideration  of  point  scattering  theory.  To  explain  the  fluctuations  of  seafloor 
reverberation,  Crowther  [1980]  modeled  the  seabed  scattering  strength  as  a  two- valued  ran¬ 
dom  variable,  and  McDaniel  [1990]  treated  seabed  reverberation  as  a  spectral  estimation 
process.  None  of  these  approaches  follow  from  a  rigorous  consideration  of  rough  surface 
scattering  theory.  An  exception  is  the  study  of  Valenzuela  and  Laing  [1970]  who  applied 
composite-roughness  theory  [Kur’ynov  1963]  to  radar  scattering  from  the  sea  surface.  They 
found,  however,  that  such  effects  were  insufficient  to  explain  the  observed  fluctuations. 

This  study  addresses  the  fluctuations  of  seafloor  reverberation  at  acoustic  frequencies  in 
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the  kHz  regime.  The  basic  assumptions  invoked  are  first  that  penetration  into  the  seabed 
with  subsequent  rescattering  may  be  neglected,  and  second  that  the  seafloor  is  slightly  rough 
so  that  scattering  may  be  treated  using  the  small  waveheight  approximation.  The  first  of 
these  assumptions  restricts  the  applicability  of  the  theory  to  the  sandy  seabeds  [Jackson 
and  Briggs  1992]  that  are  characteristic  of  the  shallow  water  areas  of  continental  shelves. 
The  second  requires  that  the  Rayleigh  parameter  be  small.  For  a  plane  wave  incident  on  a 
slightly  rough  surface  that  contains  roughnesses  of  all  length  scales,  the  small  waveheight 
approximation  predicts  that  the  backscattered  field  depends  only  on  that  surface  component 
having  a  wave  number  Kb  =  kxl‘^  where  is  the  horizontal  component  of  the  incident  wave 
vector.  However,  for  the  directional  sources  used  in  reverberation  measurements,  backscatter 
is  governed  by  a  range  of  surface  wave  numbers  about  the  Bragg  wave  number  Kb-  The  range 
of  rough  surface  wave  numbers  contributing  to  the  scatter  is  found  to  play  an  important  role 
in  determining  fluctuation  statistics. 

With  these  two  assumptions,  a  rigorous  theory  is  developed  to  predict  the  scintillation 
index  and  scattered  intensity  covariance  of  backscatter.  A  general  treatment  is  developed  for 
statistically  nonstationary  surfaces:  surfaces  for  which  the  correlation  of  the  rough  surface 
excursion  at  two  separated  points  pi  and  p2  depends  on  both  pi  and  p2  as  opposed  to  the 
stationary  case  where  the  correlation  depends  only  on  the  separation  {px  —  p2)-  The  key  to 
this  development  is  the  application  of  two  ensemble  averages:  the  first  to  that  set  of  surfaces 
having  a  fixed  or  deterministic  nonstationary  character,  and  the  second  over  all  possible  real¬ 
izations  of  nonstationary  characteristics.  The  general  theoretical  formulation  of  the  problem 
and  definition  of  the  statistical  averaging  procedures  for  stationary  and  nonstationary  surface 
repesentations  are  presented  in  section  2. 

Section  3.1  addresses  frequency  domain  scattering  from  statistically  stationary  surfaces. 
Assuming  a  simple  form  for  the  source  directivity,  the  scintillation  index  is  found  to  be 
identical  to  the  Rayleigh  value  of  unity  provided  that  the  dimensions  of  the  ensonified  area 
are  large  on  a  wavelength  scale  and  the  rough  surface  wave-number  spectrum  (the  spatial 
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Fourier  decomposition  of  the  surface  correlation  function)  is  slowly  varying  in  the  vicinity  of 
the  Bragg  wave  number.  These  two  provisions  insure  that  only  a  small  range  of  rough  surface 
wave  numbers  contribute  to  the  scatter  and  are  furthermore  consistent  with  the  geometries 
employed  in  seafloor  reverberation  measurements  and  also  with  the  applicable  measurements 
of  seabed  roughness  wave-number  spectra  [Briggs  1989].  Section  3.2  extends  the  frequency 
domain  treatment  to  statistically  nonstationary  surfaces.  Two  very  simple  representations 
of  the  nonstationary  process  are  assumed,  both  of  which  yield  similar  results  and  predict 
the  one  clear  trend,  a  significant  increase  in  fluctuations  with  decreasing  ensonified  surface 
area,  that  is  indicated  by  the  available  statistical  data  [Crowther  1980,  Chotiros  et  al  1985, 
Boehme  and  Chotiros  1988,  and  Gensane  1989]. 

In  section  4,  the  theory  is  extended  to  the  time  domain,  where  it  is  found  that  to  obtain 
results  independent  of  the  specific  form  of  the  surface  correlation  function,  the  surface  wave- 
number  spectrum  must  again  be  slowly  varying  in  the  vicinity  of  the  Bragg  wave  number.  To 
obtain  a  formulation  amenable  to  numerical  analysis,  an  additional  approximation  is  made 
that  limits  the  validity  of  the  results  to  waveforms  of  moderate  bandwidth.  Results  for  the 
scintillation  index  and  time-lagged  intensity  covariance  are  obtained  for  pure  tone  pulses, 
and  it  is  demonstrated  that  the  intensity  covariance  may  be  inverted  to  estimate  the  spatial 
structure  of  nonstationary  surface  features. 

Section  4.3  addresses  frequency  diversity  effects,  whereby  significant  reductions  of  the 
scintillation  index  have  been  experimentally  obtained  [Ulaby,  et  al  1988,  George  and  Jain 
1972]  by  increasing  the  time-bandwidth  product  of  the  waveform  used  above  its  nominal  value 
of  unity  for  a  pure  tone  pulse.  For  the  backscatter  process  considered  here,  a  frequency  di¬ 
versity  effect  is  not  predicted  for  either  stationary  or  nonstationary  surface  statistics,  given 
the  conditions  under  which  the  theory  is  valid.  This  is  an  important  finding  because  it  rep¬ 
resents  a  means  of  experimentally  differentiating  between  Bragg  scatter  and  other  scattering 
mechanisms,  such  as  volume  scatter,  that  exhibit  a  frequency  diversity  effect. 

In  section  5,  predicted  scintillation  indices  are  compared  with  the  data  of  Gensane  [1989] 
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using  a  simple  two  parameter  model  for  the  nonstationary  surface  features.  Agreement  is 
obtained  with  the  data  by  assuming  that  the  normalized  variance  of  the  surface  roughness 
wave-number  spectrum  is  approximately  unity,  and  that  the  characteristic  length  scale  of 
nonstationary  surface  features  is  on  the  order  of  1  m. 

In  section  6,  the  major  findings  of  this  study  are  summarized:  that  nonstationary  surface 
statistics  must  be  assumed  to  account  for  experimental  observations,  and  that  frequency 
diversity  effects  do  not  occur  for  Bragg  backscatter.  Other  nonstationary  features  of  the 
seafloor  that  may  affect  fluctuation  statistics  are  also  discussed. 

2.  Theoretical  development 
2.1.  General  formulation 

The  scattering  geometry  considered  is  depicted  in  figure  1,  in  which  the  acoustic  field 
incident  from  a  directional  source  is  backscattered  from  a  rough  surface.  The  backscattered 
acoustic  pressure  p{t)  obeys 

/OO 

h{T)e{t  -  T)dT  (1) 

-OO 

where  h(r)  is  the  impulse  response  of  the  scattering  channel,  and  e{t)  is  the  transmitted 
waveform.  The  impulse  response  is  related  to  the  frequency  response  H{u})  by  a  Fourier 
transform 

/OO 

dTh{T)exp(—iu)T)  (2) 

-OO 

For  a  backscattering  geometry,  the  Fraunhofer  phase  is  large,  so  that  Fresnel  phase  correc¬ 
tions  [Melton  and  Horton  1970]  may  be  neglected,  and  the  applicable  expression  for  a  slightly 
rough  surface  is 

H{uj)  =  F  j  dpiCi  exp[ia;(rj  -  axi)]  Cli  (3) 

where 

F  =  sin^  (4) 

a  =  2cos0/c 
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and  x(^)  is  a  modified  reflection  coefficient  [Jackson,  et  al,  1986],  fii  =  is  the  surface 

ensonification  function,  r,  =  2i?/c  is  the  two  way  travel  time  from  the  source  to  the  origin 
of  the  coordinate  system,  with  c  the  sound  speed  in  the  host  medium,  and  a;  the  radian 
frequency.  The  grazing  angle  6,  slant  range  R,  and  excursion  of  the  rough  surface  from  its 
mean  value  (i  =  C(Pi)  that  appear  in  these  expressions  are  depicted  in  figure  1.  In  obtaining 
this  expression  for  H{lo)  it  has  been  assumed  that  R,  9  and  x  may  be  considered  constant 
over  the  ensonified  area. 

In  the  expression  for  H{u>)  given  above,  only  its  dependence  on  the  stochastic  parameter 
Cl  has  been  made  explicit.  The  body  of  this  article  addresses  the  dependence  of  intensity 
moments  on  this  parameter.  However,  the  modified  reflection  coefficient  x(^)  may  also 
vary  with  position  through  its  dependence  on  the  local  values  of  sediment  properties.  The 
derivation  for  this  case  closely  follows  that  for  stochastic  variations  of  (  and  will  be  discussed 
later.  From  the  viewpoint  of  composite-roughness  theory,  H{u>)  also  depends  on  the  local 
large-scale  surface  slope,  through  which  the  grazing  angle  6  may  be  considered  a  stochastic 
variable.  Valenzuela  and  Laing  [1970]  investigated  this  effect  for  the  scattering  of  radar  waves 
from  the  sea  surface  and  concluded  that  it  was  insufficient  to  explain  the  experimentally 
observed  fluctuations.  Because  the  slopes  of  horizontally  large-scale  seafloor  elevations  are 
considerably  lower  than  those  of  the  sea  surface,  this  effect  is  not  considered. 

2.2.  Statistical  averages 

Our  interest  is  in  determining  the  expected  value  of  moments  of  the  received  pressure 
field,  namely  <  p{t)p*{t)  >  and  <  p{t)p*{t)p{t')p*{t')  >.  From  expressions  (1)  and  (3),  it  is 
apparent  that  these  expected  values  depend  on  the  second  and  fourth  moments  of  the  rough 
surface  excursion  <  (1(2  >  and  <  C1C2C3C4  >•  For  a  statistically  stationary  surface 

<  C1C2  >=  C(ft  -  ft)  =  J  dKW(K)  exp[i/?  ■  (p,  -  p,)\  (5) 

In  this  case  the  correlation  function  C  depends  only  on  the  separation  p\  —  p2,  with  the  surface 
roughness  wave-number  spectrum  W{K)  being  defined  as  the  spatial  Fourier  transform  of 

C. 
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If  the  rough  surface  is  statistically  nonstationary  <  C1C2  >  takes  the  form  [Papoulis,  1965] 

<  C1C2  >=  <^(^1)^2) 

With  the  change  of  variables  =  Pi-p2  and  p+  =  {pi  +  p2)l2,  a  position  dependent  surface 
roughness  wave-number  spectrum  W{K,p+)  can  be  defined 

<  C1C2  >=  C(p-,p+)  =  j  dKW{KJ+)expiK  ■  (6) 

in  analogy  to  the  definition  of  the  wave-number  spectrum  for  a  stationary  rough  surface. 
The  wave-number  spectrum  W{K)  averaged  over  all  nonstationary  realizations  of  the  rough 
surface,  where  the  overbar  is  used  to  designate  this  average,  is  defined  as 

_  ry+  dp^  ,  -  ^  . 

Let  us  now  return  to  the  problem  at  hand,  the  estimation  of  such  quantities  as 

<  >=  W?  j  dp-dp^dKW{k,pj^)9.i9,2exp{i[K  ■  p-  -  ua{xi  -2:2)]} 

The  ensonification  factors  fli  and  fl2  in  this  expression  determine  the  extent  of  the  integra¬ 
tion  over  This  integration  then  yields  a  local  average  of  the  position  dependent  wave- 
number  spectrum  over  the  ensonified  area.  Thus,  to  obtain  the  expectation  of  H{u)H*{ui) 
for  all  nonstationary  representations  of  the  rough  surface,  an  additional  ensemble  average 
corresponding  to  that  performed  to  obtain  W{K)  is  necessary 

<  H{u)H*{u:)  >=  \F\^  J  dp.dp+dkW{k)nin2exp{i[k  •  p.  -  uja{xi  -  2:2)]} 

The  procedure  described  above  for  determining  the  expected  moments  of  the  rough  sur¬ 
face  excursion  entails  two  ensemble  averages.  The  first  is  over  all  realizations  of  the  surface 
having  the  same  nonstationary  characteristics.  An  example  of  such  a  set  of  surfaces  is  an 
ensemble  in  which  the  rms  roughness  varies  deterministically  with  position.  The  second 
ensemble  average  is  over  all  nonstationary  realizations  of  the  suface.  For  the  example  con¬ 
sidered,  this  second  average  is  over  all  possible  variations  of  the  rms  roughness  with  position. 
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Obtaining  second  monaents  of  the  scattered  intensity  entails  evaluating  the  fourth  moment 
of  the  rough  surface  excursion  <  C1C2C3C4  >•  If  C  is  a  Gaussian  random  variable 

<  C1C2C3C4  >=<  C1C2  X  C3C4  >  +  <  C1C3  ><  C2C4  >  +  <  C1C4  X  C3C2  >  (7) 

It  follows  from  the  result  obtained  above  for  <  C1C2  >  that  the  ensemble  average  over  all 
nonstationary  surfaces  of  a  representative  term  on  the  right  hand  side  of  (7)  takes  the  form 

<  C1C2  X  C3C4  >=  / dKdK'w(K, ■  {pi  -  P2)  +  iii'  ■  (P3  -  Ml 

where  =  (^3  +  P4)/2. 

If  W(K)  for  a  nonstationary  surface  is  chosen  equal  to  W(K)  for  a  stationary  surface, 
then  the  mean  intensity  scattered  in  the  two  cases  will  be  identical.  However,  the  above 
expression  for  the  fourth  moment  of  (  depends  on  the  correlation  function  of  W(K,^  at 
disjoint  positions,  so  that  the  second  moment  of  the  intensity  will,  in  general,  differ  in  the 
two  cases. 

3.  Frequency  domain  scattering 

It  is  clear  that  obtaining  numerical  results  from  the  theoretical  formulation  of  section 
2  presents  some  difficulty.  It  is  hence  useful  to  first  examine  scattering  in  the  frequency 
domain  to  establish  the  validity  of  the  various  approximations  that  are  needed  to  obtain 
results  amenable  to  numerical  analysis. 

S.l.  Statistically  stationary  surfaces 

To  obtain  results  for  the  statistically  stationary  case,  it  will  be  assumed  that  the  projector 
ensonification  function  ft  is  Gaussian 

fli  =  exp[-x^/(Ax)2  -  yll{Ay)^]  (8) 

where 

Ax  =  gRta.n{(f>xl2)l  sinO 
Ay  =  gRta,n{(j)yf2) 
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with  <f>x  and  (f)y  the  projector  half-power  beam  widths  along  the  x  and  y  axes,  respectively, 
and  g  =  [2/ln(2)]^/^.  The  ensonified  area  in  this  case  is  given  by  7rAxAy/2. 

For  this  form  of  the  ensonification  function,  <  /(w)  >=<  >  is  given  by 

<  I{u)  >=  |Fp  J  dp+dp.dKW{K)  exp[i{K  ■  p.  -  c^ax.)] 

2x\  +  x'il2  2yl  +  yll2 
- (A^)^ - 

where  p-  =  (x_,j/_)  and  p+  =  {x+,y+).  The  integrals  over  and  p+  in  this  expression 
may  be  immediately  performed  to  obtain 

<  /(cu)  >=  \F\V{AxAyYj  dKW{K)exp[-{K,  -  au:)\Axfl2  -  Kl{Ayfl2] 

The  remaining  integral  over  K  is  an  average  of  W{K)  about  the  Bragg  wave  number  Kb  = 
(aa;,  0).  With  the  definition 

Wb  =  AxAy  J  exp[-(/4  -  c^f{Axfl2  -  K^,{Ayf/2] 

the  mean  scattered  intensity  takes  the  form 

<  I{uj)  >=  2Tr^\FY{AxAy)WB  (9) 

It  should  be  noted  that  for  wave- number  spectra  that  are  strongly  varying  functions  of  K,  Wb 
may  differ  significantly  from  W{Kb)-  With  the  same  approximations,  the  second  moment 
of  the  intensity  is  given  by 

<  P{u;)  >=  2  <  I{u)  (10) 

+{\FYTv\AxAyY  J  dKW{K)exp[-{Kl  +  aV)(Ax)V2  -  Kl{AyY/2]Y 

The  expressions  developed  above  for  the  first  and  second  moments  of  the  backscattered 
intensity  will  be  evaluated  for  two  assumed  forms  of  the  wave- number  spectrum  W{K).  The 
first  is  the  Gaussian  form  that  is  often  assumed  in  scattering  problems  because  of  its  ease  of 
numerical  implementation 

W(K)  =  Sexfl-(Kl  +  Kl)L^] 
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where  Z  is  a  characteristic  length,  and  S  a  scale  factor.  For  this  wave-number  spectrum, 

Wb  =  S{[l+2{L/Axy][l  +  2{L/Ayy]}-'^/^exp{-{aujL)y[l+2{L/Axy]} 

It  is  evident  from  this  expression,  that  for  L  >  Ax,  Wb  will  greatly  exceed  W{Kb)  = 
S exp[-(ac<;Z)^].  In  such  cases,  a  broad  range  of  wave  numbers  contributes  to  Wb,  and  the 
scattering  process  differs  from  classical  view  of  Bragg  diffraction  as  wave-number  selective. 

For  the  Gaussian  spectrum,  the  scintillation  index  (7^  =  [<  >  /  <  I{UJ  >2]  -  1 

takes  the  form 

CT^  =  1  +  exp{-(aa)Ax)^/[l  +  2(Z/Ax)^]} 

It  is  clear  from  this  result  that  for  a  Gaussian  spectrum,  cr^  depends  strongly  on  frequency, 
a  dependence  that  has  not  been  observed  experimentally.  This  result  is  also  independent  of 
Ay  the  cross-range  extent  pf  the  ensonified  area. 

The  second  surface  wave-number  spectrum  that  will  be  considered  is  more  typical  of 
those  associated  with  the  seafloor  [Briggs  1989,  Fox  and  Hayes  1985,  Goff  and  Jordan  1988] 

W{K)  =  S/[l  +  L\Kl  +  K^^)]  (11) 

The  numerical  evaluation  of  Wb  and  for  this  wave-number  spectrum  may  be  simplified 
by  the  use  of  the  formula  [Gradshteyn  and  Ryzhik  1980] 

r®®  duexp[— (nu)^] 

Jo  " -(1,2  +  ^2)  =  exp[(/?^)  ]/(2/?) 

where  $  is  the  probability  function,  to  perform  the  integrations  over  Ky. 

The  scintillation  index  predicted  for  the  power-law  spectrum  is  compared  with  that  for 
the  Gaussian  wave-number  spectrum  in  figure  2.  In  this  example,  the  grazing  angle  0  is  30°, 
the  frequency  10  kHz,  and  the  characteristic  length  L  is  20  m.  For  the  narrow  beam  patterns 
assumed  (j)^  —  <f>y  =  1°,  the  dimensions  of  the  ensonified  area  do  not  exceed  L  even  for  the 
longest  range  considered.  The  results  presented  in  figure  2  parallel  those  obtained  by  Yang, 
Fennemore  and  McDaniel  [1992]  for  forward  scattering.  They  found  that  strong  fluctuations 
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occured  only  when  the  source  or  receiver  was  in  the  near  field  of  the  surface,  and  that  these 
fluctuations  were  stronger  for  a  Gaussian  wave-number  spectrum  than  for  a  quasipower-law 
spectrum.  They  also  found  that  the  introduction  of  an  ensonification  function  led  to  higher 
fluctuations. 

For  the  extreme  example  considered  in  figure  2,  Wb  predicted  by  the  Gaussian  wave- 
number  spectrum  exceeds  W{Kb)  by  many  orders  of  magnitude.  The  ratio  WbIW{Kb)  is 
shown  for  the  power-law  spectrum  in  figure  3.  Only  at  the  very  shortest  ranges  does  this 
ratio  differ  significantly  from  unity.  It  is  apparent  from  the  results  of  this  subsection  that  if 
we  restrict  our  interest  to  wave-number  spectra  that  are  slowly  varying  in  the  neighborhood 
of  the  Bragg  wave  number,  and  to  ensonification  functions  such  that  Aa;  and  Ay  are  large 
on  a  wavelength  scale,  the  scintillation  index  will  be  unity  for  scattering  from  stationary 
surfaces. 

3.2.  Nonstationary  surfaces 

For  the  nonstationary  case,  the  derivation  closely  follows  that  for  the  stationary  case. 
With  the  introduction  of  a  Gaussian  beam  pattern,  the  mean  intensity  for  a  nonstationary 
representation  of  the  rough  surface  may  readily  be  reduced  to  the  form 

<  I{uj)  >=  47r2|F|2|  JpVWB(pV)exp[-24/(Ax)2  -  2yll{Ayf] 

where 

Ws(p+)  =  Ax  Ay  j  ^w[-{Kx  -  OLuf{Axf  {2  -  Kl{Ayf  12] 

An  ensemble  average  over  all  nonstationary  surfaces  then  yields 

<  /(w)  >  =  2'K^\F^{AxAy)WB 

where  it  has  been  assumed  that  W b{p+)  =  Wg. 

It  is  clear  from  the  results  of  setion  3.1  that  for  cases  of  practical  interest  <  /^(w)  >  may 
be  approximated  by  2  <  I{io)  >^.  Performing  an  average  over  all  nonstationary  representa¬ 
tions  of  the  surface  then  yields 

<  P{iJL>)  >  =  2<  I{u})  ><  /(w)  > 
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<  PM  >  =  2(2^)‘|f  |‘  J  dp^dplWBlMWB(p[) 

X  exp[-2(x^  +  x^)/{Axy  -  2(y+  +  y'^)/{Ayf] 

Because  Wb(p+)  and  Wb{p[)  entail  an  average  over  a  narrow  band  of  wave  numbers  about 
the  Bragg  wave  number  for  wave-number  spectra  of  interest,  their  correlation  function  will 
be  assumed  to  depend  on  position  alone. 

To  evaluate  the  expression  for  <  P{>^)  >  the  two  position  wave-number  correlation  func¬ 
tion  is  assumed  to  be  of  the  general  form 

Wb{p+)Wb{p!,)  =  W^s[i+  'fCwip^  -  P+)]  (12) 


where  7  is  the  normalized  variance  of  the  wave-number  spectrum  and  Cw  is  a  normal¬ 
ized  correlation  function.  Using  this  form  for  the  correlation  function  in  the  expression  for 


<  P{oj)  >,  thence  making  the  change  of  variables  u  =  p+  —  p|,  U  =  {p+  +  p|)/2,  the 
scintillation  index  may  be  reduced  to 


=  1  +  / 


27  du 


u 


— - — —Cwiu)  exp[-  ^7- 
nAxAy  (Ax)^  {AyY 


y _ ■ 

..\i- 


Crowther  [1980]  modeled  scattering  from  a  seafloor  consisting  of  random  patches  taking 
on  only  two  values  of  the  scattering  strength:  a  low  scattering  strength  with  probability 
s  and  a  high  scattering  strength  with  probability  1  —  s.  Assuming  a  ratio  of  high  to  low 
scattering  strengths  z^,  he  found  7  =  s(l  —  s)(i/  —  1)^/[1  -|-  s{v  —  1)]^  and  that  Cw  was 
exponentially  distributed 

Cw{u)  =  exp{-\u\/ Lw) 

where  Lw  characterizes  the  patch  size.  Inserting  this  correlation  function  in  the  expression 
for  cr^  and  making  the  change  of  variables  =  rcos^’,  Uy  =  rsin0,  yields 


where  /(V’)  =  -ip I (Ax)^  +  cos^  ^/(Aj/)^.  The  integration  over  r  in  this  expression  may 
be  readily  performed.  The  remaining  integral  over  !/>  can  then  be  evaluated  numerically. 
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To  examine  the  sensitivity  of  the  scintillation  index  to  the  assumed  form  of  the  correlation 
function  the  Gaussian  case  will  also  be  considered 

Cw{u)  =  exp[— (u^  +  Uy)/ L^]  (13) 

For  this  correlation  function, 

<7^  =  1  +  27/{[1  +  (Ax/Ti^)^][l  + 

It  is  apparent  from  this  result  that  1  <  <  1  +  27,  with  the  maximum  value  being  obtained 

when  Lw  is  much  larger  than  the  dimensions  of  the  ensonified  area. 

The  results  obtained  above  for  the  scintillation  index  display  no  explicit  dependence 
on  frequency.  However,  such  a  dependence  is  implicit  because  Wb{p-^)  and  its  correlation 
depend  on  the  Bragg  wave  number  Kb  =  (oiWjO). 

The  scintillation  indices  predicted  for  the  exponential  and  Gaussian  correlation  functions 
are  compared  in  figure  4  which  shows  their  dependence  on  slant  range  for  a  fixed  grazing  angle 
30°  and  beam  pattern  (l>x  =  <f>y  =  10°.  In  this  example,  7  =  .25  and  Lw  =  20  m.  The  results 
for  the  exponential  and  Gaussian  correlation  functions  are  in  close  agreement  and  approach 
unity  as  the  slant  range  increases,  and  the  dimensions  of  the  ensonified  area  approach  the 
characteristic  length  Lw-  This  decrease  in  the  scintillation  index  is  more  evident  in  figure  5 
which  shows  corresponding  results  for  Lw  =  5  m.  In  this  figure,  the  results  obtained  for  the 
two  correlation  functions  are  again  in  close  agreement. 

Figure  6  shows  the  dependence  of  on  4>x  lor  a  fixed  grazing  angle  of  30°  for  two  values 
of  the  characteristic  length:  the  remaining  parameters  are  those  used  in  figure  5.  Again, 
the  predicted  scintillation  index  decreases  as  the  size  of  the  ensonified  area  increases.  In 
this  example,  clear  differences  in  the  predictions  of  the  Gaussian  and  exponential  models  are 
evident  for  very  small  beamwidths. 

4.  Scattering  in  the  time  domain 

In  most  backscattering  experiments  the  resolution  in  range  is  determined  by  the  length 
of  the  transmitted  pulse  and  the  effect  of  the  beam  pattern  on  range  resolution  is  negligible. 
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Hence,  in  obtaining  results  in  the  time  domain  ^{x,y)  will  be  replaced  by  n{y)  in  the 
analysis. 

An  additional  approximation  will  be  made  to  simplify  the  analysis:  that  the  scatter¬ 
ing  is  frequency  independent.  With  this  assumption  the  function  F  is  approximated  by 
Lolxi^)  sin^  6 1 [tt R'^ c^)  where  uq  is  the  center  frequency  of  the  transmitted  pulse.  While  this 
approximation  is  consistent  with  those  made  to  approximate  other  factors  in  F  that  do  not 
appear  in  the  argument  of  an  exponential  by  their  mean  values,  it  is  evident  that  it  limits 
the  validity  of  the  results  to  waveforms  of  moderate  bandwidth. 

Following  the  general  formulation  developed  in  section  2,  one  finds 

h{T)  =  F  j  dxxdyi8{T  -t,-  a;iQ:)exp[-j/i/(Ay)^]Ci 

and  consquently  that 

p{t)  =  F  j  dxidyiCi  exp[-yl/{Ayy]e[t  -  (r*  +  xja)] 

For  statistically  nonstationary  surfaces,  it  follows  that 

<  p{t)p-{f)  >=  \F\^  J ip.dp+dKW(K,Mexpl-{yl  +  yl)/{Ayf] 

X  exp{iK  •  p-)e[t  —  (r^  -|-  xiC()]e*[t'  —  (r^  +  X20c)]  (14) 

Irrespective  of  the  frequency  independent  approximation  used  to  obtain  (14)  the  time- lagged 
intensity  covariance  is  given  by 

<  >=<  p{t)p*{t)  X  p{t')p*{t')  >  +  <  p{t)p*{t')  ><  p{t')p*{t)  >  (15) 

where  the  only  assumption  made  to  obtain  this  result  is  that  an  additional  term  <  p{t)p{t')  > 
<  p*{t)p*(t')  >  appearing  in  this  expression  may  be  neglected,  as  demonstrated  in  section  3. 

Before  considering  specific  waveforms,  it  is  worthwhile  obtaining  general  results  for  sta¬ 
tistically  stationary  surfaces  to  determine  if  any  approximations  that  simplify  the  analysis 
may  be  obtained.  First  consider 

<  /i(r)/i*(r')  >=  |Fp  f  dpidp2S{T  —  Tg  —  Xia)8{T'  —  r,  —  X2a) 
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xexpl-(!/;  +  !,|)/(Aj,)^lC(p-i-ft,)  (16) 

Replacing  C[p\  —  p2)  by  its  Fourier  decomposition  in  wave-number  space,  the  result  obtained 


is 

<  h(T)h.-(T')  >=  miaf  j  dy+dy.dKW{K) 

X  expl-(2!,J  +  yl/2)/{^yf  +  iKAr  -  T')/a  +  tK,y.] 

Performing  the  integrals  over  y+  and  j/_  in  this  expression,  then  yields 

<  h{T)h*{T')  >=  7r{Ay\F\/afj  dKW{K)exp[iK,{T  -  r')/a  -  K^{Ay)V2] 

This  integral  may  be  readily  evaluated  numerically.  Normalized  results  are  shown  in  figure 
7  for  the  power-law  spectrum  (11)  for  two  values  of  the  characteristic  length  L.  It  is  evident 
from  this  figure,  that  the  temporal  covariance  of  the  impulse  response  is  a  slowly  varying 
function  of  lag  time  for  the  time  scales  of  interest.  It  is  also  apparent  that  this  function  is 
dependent  on  the  characteristic  length  L,  as  is  strongly  suggested  by  the  representation  (16). 

Because  consideration  of  the  time-lagged  covariance  of  the  impulse  response  fails  to  yield 
any  obvious  simplifications,  we  next  address 

<  p{ts  +  Ai/2)p*(rj  -  At/2)  >=  \F\^  J  dp-dp+dKW{K)exp{iK  ■  p_) 

X  exp[-(2y2  +  yl/2)/iAy)M^t/2  -  x,a)e*i-At/2  -  x^a) 

where  t  and  t'  have  been  replaced  by  -|-  A^/2  and  -  At  12,  respectively,  and  stationary 
surface  statistics  have  been  assumed.  Performing  the  integrals  over  j/+  and  j/_,  and  substi¬ 
tuting  e{At/2  —  x\a)  =  e(Af/2  —  xio;)  exp[iu;o(At/2  —  xia)],  where  e  is  the  complex  signal 
envelope,  along  with  the  corresponding  relation  for  e{—At/2  —  X2(x),  into  this  expression 
yields 

<  p{ts  -1-  At/2)p*{Ts  -  At  12)  >=  7r|Fp(Aj/)^  J  dUdKydxidx2W{U  +  aujo,  Ky) 

X  exp[—Ky{Ay)^/2  +  iU{xi  —  X2)  +  iu}oAt]e{At /2  —  Xia)e*{—At/2  —  X2a)  (17) 
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where  U  =  ^  it  is  now  assumed  that  W{U  +  au>o,Ky)  varies  so  slowly  with  U 

that  this  dependence  may  be  neglected,  the  integral  over  U  yields  27r^(xi  —  X2),  and 

<  p{ts  +  At/2)p*(Ts  -  At/2)  >=  27r^lF|^(A2/)^exp(zu;oAt) 

X  J  dx+e{Atl2  -  x+a)e*{-At/2  -  x+a)  j  dKyW{a(jjQ,Ky)eY.p[-Kl{Ayf  12]  (18) 

The  analagous  approximation  for  nonstationary  surfaces  takes  the  form 
<  p(r^  +  Atl2)p*{T,  -  At/2)  >=  (27r)2|T’|^exp(ia;oAt) 

X  j  dx+e{At/2  —  x+a)e*{—At/2  —  x^a)  j  dy^WB{p+)exp[—2y^/{Ay)^]  (19) 

where  Wb{p+)  is  defined  as 

Wb{M  =  Ayl{2-Kfl'^  j  dKyW[{ocuo.Ky),p+\txp-Kl{Ayf  12 

It  is  clear  that  (18)  could  also  be  obtained  by  assuming  that  <  h{T)h*{T')  >  of  (16)  contains 
a  factor  8{t  —  t')  which  is  the  case  for  point  scatterers  [Ol’shevskii  1967].  However,  we  have 
seen  that  <  h{r)h*{T')  >  is  in  general  a  slowly  varying  function  of  r  —  r'  and  hence  such 
an  assumption  is  useful  for  rough  surface  scatter  only  if  the  surface  roughness  wave-number 
spectrum  is  a  slowly  varying  function  of  Kx  near  Kx  =  ohhq. 

4.1.  Pure  tone  pulses 

The  expressions  above  will  first  be  evaluated  for  a  Gaussian  weighted  pure  tone  pulse 

e{t  -  t)  =  exp[iwo(i  -  t)  -  {t  -  t)^/T‘^]  (20) 

for  which  the  half-power,  or  effective,  pulse  length  Tp  is  given  by  (21n2)^'^^r.  First,  we 
consider  the  scintillation  index,  which  entails  evaluation  of  <  p{ts)p*{ts)  >.  Proceeding 
directly  to  the  case  of  statistically  nonstationary  surfaces  one  readily  finds 

<  P(r,)p-(r,)  >=  teVP  /  dp+WaiM  expl-2yl/{Ayy  -  x\cfl1^\ 

where 

WBih)  =  {AyT/2na)  J  dKW{K,  p+)exp[-Kl{Ay)y2  -  (Kx  -  {2a^)] 
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On  comparing  this  result  with  the  comparable  result  for  <  /(w)  >  obtained  in  section  3.2  for 
scattering  of  continuous  wave  signals  from  nonstationary  surfaces,  it  is  apparent  that  with 
the  substitution  T/a  =  Ax,  the  two  expressions  are  identical.  It  thus  follows  that  the  results 
presented  in  figures  4  -  6  for  the  scintillation  index  are  applicable  also  to  the  backscatter  of 
Gaussian  weighted  pure  tone  pulses. 

To  compute  the  time-lagged  intensity  covariance,  expressions  are  needed  for  <  p{t)p*{t)  > 
and  <  p{t')p*{t')  >  or  for  <  /(r^  i  At/2)  >.  For  the  pure  tone  pulse  (20) 

<I(r,tAt/2)>={2,r\F\f  Jdp^%(p+)  (21) 

X  exp[-24/(A!;)*  -  2xla^/T'‘  t  2Ate+Q:/r’  -  (Atf/{2T^)] 

Taking  an  ensemble  average  over  all  nonstationary  surface  representations 

<  /(r,  i  At/2)  >=  27r^\F\^Ay{T/a)WB 

yields  a  result  that  is  independent  of  time,  where  Wb  =  Wb{M- 

Because  the  results  presented  in  section  3  display  only  a  weak  dependence  on  the  form  of 
the  correlation  function  assumed  for  Wb{p+)Wb{p!1)  in  (12),  a  Gaussian  correlation  function 
(13)  will  be  assumed.  For  the  Gaussian  case,  the  first  term  contributing  to  the  time-lagged 
intensity  covariance  (15)  takes  the  form 

<  I{t,  + At/2)  At/2)  >=  {2Tr^\F\^WBAyT/a)^ 

x{l +'f'exp[-{At/aLw)^/{'i- +T'^/o/^Lly)]}  (22) 

where  7'  is  defined  as 

y  =  7[1  +  (AylLwfV'l^l  +  (TlaLwf]-'!^ 

For  a  stationary  surface,  this  contribution  to  the  intensity  covariance  is  independent  of  lag 
time.  In  the  case  of  nonstationary  surfaces,  a  time  dependence  is  predicted  that  depends 
strongly  on  the  characteristic  patch  size  Lw  A  result  of  this  general  form  was  first  obtained 
by  Crowther  [1980]. 
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Using  approximation  (19),  the  second  contribution  to  the  time-lagged  intensity  covariance 
is 

1  <  p(r^  -H  Af/2)p*(rs  -  Ai/2)  >  p  =  {2'K^\F\^WB^yT j af{\  +  7')exp[-(A^)VT’^]  (23) 

The  dependence  of  this  contribution  to  the  intensity  covariance  on  lag  time  is  the  same 
for  both  stationary  and  nonstationary  surface  statistics,  the  only  difference  in  the  two  cases 
being  the  scale  factor  1  +7'.  On  comparing  (22)  and  (23),  it  is  evident  that  the  nonstationary 
contribution  of  (22)  decays  more  slowly  with  lag  time  than  that  of  (23)  and  hence  provides 
a  means  of  estimating  Lw- 

Assuming  that  the  wave-number  spectrum  in  (22)  is  a  slowly  varying  function  of  Wb 
in  this  expression  may  be  replaced  by  Wb-  We  then  consider  the  normalized  time- lagged 
intensity  covariance  C{At)  where 

C{At)  =<  /(r.-H  Ai/2)/(r,-AV2)  >  /[<  /(r,)  >]2 

In  figure  8,  the  normalized  intensity  covariance  is  shown  as  a  function  of  lag  time  for  station¬ 
ary  and  nonstationary  surfaces.  This  figure  presents  predictions  for  a  slant  range  of  50  m, 
grazing  angle  of  30°,  and  horizontal  beamwidth  of  10°.  Results  are  shown  for  a  pulse  length 
of  1  msec,  for  three  characteristic  patch  sizes.  For  the  value  of  7  =  .5  and  the  scattering 
geometry  assumed  in  this  figure,  only  for  Lw  =  5  m,  will  the  slowly  varying  time  dependent 
factor  in  (22)  be  clearly  apparent,  given  that  the  validity  of  these  results  requires  At  Tg. 

Figure  9  shows  the  normalized  intensity  covariance  predicted  for  a  scattering  geometry 
of  larger  scale  in  which  a  slant  range  of  500  m  and  pulse  length  of  10  msec  are  assumed.  A 
horizontal  beamwidth  and  grazing  angle  of  10°  were  used  to  obtain  results  for  a  stationary 
surface  7  =  0  and  a  nonstationary  surface  with  7  =  1.  For  this  pulse  length,  nonstationary 
features  having  a  charateristic  length  Lw  <  5  m  are  not  experimentally  resolvable.  However, 
the  decay  of  the  covariance  due  to  large  scale  nonstationary  features  permits  an  estimation 
of  Lw  in  such  cases. 
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4-2.  Spectral  representation  of  nonstationary  processes 

While  the  examples  presented  in  figures  9  and  10  provide  some  insight  into  the  depen¬ 
dence  of  the  intensity  covariance  on  experimental  and  seafloor  parameters,  they  entail  the 
assumption  that  Cw  depends  on  a  single  length  scale  Lw  A  more  general  form  for  this 
correlation  function  is 

Cw{^  =  j  dKS{K)  expiK  •  p 

For  a  pure  tone  pulse,  the  contribution  to  the  intensity  covariance  of  (22)  then  takes  the 
form 

<  /(r, -H  At/2)  X /(r,  -  At/2)  >=  {2ir^\F\^WBAyT/aY 

x{l  -f  (7/47r2)  J  dKSiK)  exp[-K^^{Ay)yA  -  +  iAT.At/a]}  (24) 

With  the  change  of  variables  At  =  2A/c,  (24)  may  be  Fourier  transformed  with  respect  to 
A  via  T{k^)  =  / dkexp{-iK^K)f{K).  Assuming  a  record  of  sufficient  length,  there  results 

:F(/c.)  =  {jac/iir){2n^\F\'‘WBAyT/a)'‘  j dK,S{acKj2,K,) 

X  exp[-c^rV^/4  -  Kl(AyYH]  (25) 

This  procedure  thus  provides  a  damped  and  averaged  estimate  of  the  spatial  decomposition  of 
Cw-  For  (25)  to  provide  a  reliable  estimate  of  5( Ala;,  0),  it  is  required  that  Ay  >  Amax  where 
Amax  is  the  maximum  length  scale  of  interest  so  that  averaging  effects  will  be  unimportant. 
In  addition,  for  the  rolloff  with  increasing  in  (25)  to  be  negligible,  we  require  T  <  Amin/ ttc 
where  Amin  is  the  minimum  length  scale  of  interest. 

Chotiros,  et  al  [1985]  Fourier  transformed  detrended  intensity  records  and  obtained  re¬ 
sults  that  differ  from  predictions  [McDaniel  1990]  for  1  m  <  A  <  10  m.  Their  data  were 
acquired  using  a  pulse  length  of  1  msec  and  Ay  «  3.5  m,  so  that  neither  of  the  criteria  above 
are  fulfilled,  and  the  data  cannot  be  inverted  to  estimate  5'(/C,0). 
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4-3.  Frequency  diversity 

In  both  microwave  and  optical  scattering  from  rough  surfaces  a  frequency  diversity  or 
frequency  averaging  effect  has  been  experimentally  observed  [Ulaby  et  al  1988,  George  and 
Jain  1972].  For  backscatter  this  effect  can  be  characterized  by  the  parameter  N  «  2D/S,f  jc 
[Ulaby  et  al]  where  D  is  the  slant-range  resolution  of  the  detection  system,  and  A/  is  the 
bandwidth  of  the  transmitted  waveform.  Experiments  employing  broadband  signals  have 
been  used  to  verify  that  cr^  w  l/iV  for  scattering  from  sufaces  that  are  homogeneous  on  the 
scale  of  the  sensor’s  footprint,  so  that  as  the  signal  bandwidth  is  increased  from  its  nominal 
value  for  a  pure  tone  pulse  A/  =  l/T,  the  scintillation  index  decreases.  It  is  of  interest  to 
determine  if  such  an  effect  is  predicted  in  the  case  of  Bragg  backscatter. 

First  let  us  consider  scattering  from  statistically  stationary  surfaces.  If  (15)  is  valid  then 

<  P{t)  >=  2  <  lit) 

so  that  (T^  =  1  regardless  of  the  transmitted  waveform  and  no  frequency  diversity  effect 
is  predicted.  The  validity  of  (15)  requires  that  the  dimensions  of  the  scattering  area  be 
large  on  the  scale  of  the  wavelength  of  the  ensonifying  signal  and  that  the  surface  roughness 
wave-number  spectrum  vary  slowly  in  the  vicinity  of  the  Bragg  wave  number  so  that  the 
scattering  process  is  highly  wave-number  selective.  Recall  also  that  (15)  is  valid  irrespective 
of  the  frequency  independent  assumption. 

For  scattering  from  nonstationary  rough  surfaces,  we  invoke  approximation  (19)  and 
assume  a  waveform 

e(f)  =  exp[i(c<;oi  +  bF)  —  t^  jT"^] 

It  may  be  readily  verified  that  in  this  case  the  result  obtained  is  independent  of  the  choice  of 
the  parameter  b  so  that  the  results  reduce  to  those  obtained  for  the  pure  tone  pulse  of  (20)  and 
again  a  frequency  diversity  effect  is  not  predicted.  In  this  case,  two  additional  assumptions 
have  been  made:  that  the  scatter  is  frequency  independent,  and  that  the  approximation  of 
(17)  by  (18)  is  valid. 
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It  is  appropriate  at  this  point  to  examine  in  detail  the  validity  of  approximation  (18).  For 
a  pure  tone  pulse,  (17)  and  (18)  may  be  integrated  assuming  wave-number  spectra  of  Gaus¬ 
sian  and  power- law  form.  In  the  Gaussian  case  (17)  yields  a  dependence  on  At  of  the  form 
exp[— (At)^/2(r^  -|-2a^L^)]  while  approximation  (18)  yields  a  dependence  exp[— (Af)^/2T^]. 
It  is  clear  that  only  if  T  l:|>  aL  is  (18)  valid  for  this  assumed  form  of  the  wave-number 
spectrum.  Results  are  shown  for  the  power-law  spectrum  of  (11)  in  figure  10,  for  a  slant 
range  of  50  m,  grazing  angle  of  30°,  horizontal  beamwidth  of  10°  and  frequency  of  10  kHz. 
The  three  curves  shown  are  the  normalized  level  predicted  by  (18),  the  relative  magnitude  of 
the  predictions  of  (18)  to  those  of  (17),  and  the  difference  in  phase  between  the  predictions 
of  (18)  and  (17).  It  is  evident  that  for  the  power-law  wave-number  spectrum,  approximation 
(18),  and  hence  (19),  is  valid  even  for  the  values  of  T  =  1  msec  and  L  =  100  m,  used  in  this 
example,  for  which  T  <  aL.  Thus,  within  the  additional  limitation  imposed  by  the  frequency 
independent  assumption,  no  diversity  effect  is  predicted  for  scatter  from  nonstationary  rough 
surfaces. 

The  reason  for  the  lack  of  a  frequency  diversity  effect  for  the  Bragg  backscattering  pro¬ 
cesses  addressed  in  this  study  is  not  intuitively  evident.  It  is  thus  worthwhile  seeking  a 
simple  physical  explanation  for  this  finding.  In  the  case  of  point  scatterers,  the  diversity 
effect  occurs  because  each  irradiated  scatterer  returns  a  replica  of  the  transmitted  pulse.  For 
Bragg  scatter,  however,  the  scatter  at  frequency  wi  is  governed  by  surface  features  having 
a  wave  number  =  auji,  and  at  frequency  u}2  by  those  having  a  wave  number  =  aL02‘ 
Thus,  in  the  case  of  Bragg  backscatter,  frequency  diverse  waveforms  increase  the  effective 
number  of  “scatterers”,  which  has  no  effect  since  this  number  is  already  large.  For  the 
same  reason,  the  effects  of  frequency  diversity  may  not  be  fully  realized  for  other  resonant 
scattering  processes,  such  as  scatter  from  the  resonant  microbubbles  that  are  responsible  for 
high-frequency  sea  surface  backscatter. 

Frequency  averaging  effects  have  been  experimentally  observed  for  both  volume  scatter 
and  backscatter  from  rough  surfaces  at  elevated  grazing  angles.  Thus,  the  absence  of  a 
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frequency  diversity  effect  is  an  important  feature  of  Bragg  scattering  processes  permitting 
their  clear  identification.  It  is  noteworthy  that  the  results  of  two  experiments  reported  by 
Chotiros,  et  al  [1985]  and  Boehme  and  Chotiros  [1988],  where  the  parameter  iV  took  on  values 
of  one  and  four  at  both  sites  where  measurements  were  conducted,  display  no  dependence 
of  the  scattered  field  statistics  on  the  bandwidth  of  the  transmitted  waveform. 

5.  Comparison  with  experiment 

The  data  with  which  the  theoretical  predictions  will  be  compared  were  acquired  by  Gen- 
sane  [1989]  at  two  sites  in  an  area  characterized  by  a  sandy  seafloor.  At  one  of  the  sites, 
the  seafloor  was  described  as  flat,  while  at  the  second,  sand  ripples  with  a  spatial  period  of 
1.2  m  were  observed.  A  parametric  source  having  a  nominal  horizontal  beamwidth  of  3°  at 
the  frequencies  of  interest  was  used  in  these  measurements  which  were  conducted  along  four 
bearings  at  each  site.  The  height  of  the  projector  above  the  seabed  was  adjusted  so  that  the 
slant  range  to  the  ensonified  area  varied  between  20  and  40  m.  The  projected  signals  were 
selected  to  obtain  data  for  a  wide  range  of  frequencies,  pulse  lengths,  and  signal  bandwidths. 

Our  interest  is  in  the  pulse  length  dependence  of  the  scintillation  index  reported  for 
gated  pure  tone  pulses  at  frequencies  of  20  and  40  kHz  over  the  grazing  angle  range  of  4  to 
20°.  These  indices  were  determined  by  first  analyzing  the  data  to  determine  the  scattering 
strength  and  then  detrending  the  data  to  remove  the  general  dependence  of  the  scattering 
strength  on  grazing  angle,  a  procedure  described  in  detail  by  Chotiros,  et  al  [1985].  Further 
analysis  was  then  performed  to  determine  the  probability  distribution  and  various  moments 
of  the  detrended  data. 

The  measured  scintillation  indices  shown  in  figure  11  represent  an  average  over  the  data 
at  the  two  sites.  The  data  for  pulse  lengths  of  .25  and  1  msec  were  acquired  at  a  frequency 
of  40  kHz.  The  20  kHz  data  plotted  for  a  pulse  length  of  1.5  msec  is  the  reported  average  of 
scintillation  indices  measured  using  1  and  2  msec  pulses.  The  theoretical  curves  presented 
for  comparison  assume  Cw  of  (12)  is  Gaussian,  in  which  case 

=  1  +  27[1  +  +  (Ay/W)^]-!/^ 
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for  a  Gaussian  weighted  pure  tone  pulse,  and  it  is  assumed  that  this  expression  is  applicable 
for  a  gated  pure  tone  of  duration  equal  to  the  effective  pulse  length  of  the  Gaussian  weighted 
pulse.  An  average  slant  range  of  30  m  and  grazing  angle  of  12°  were  used  in  evaluating  cr^. 

As  is  evident  with  the  logarithmic  time  scale  employed  in  figure  11,  cr^  exhibits  a  “knee” 
whose  position  is  strongly  dependent  on  the  characteristic  patch  size  Lw  assumed,  and  to 
obtain  agreement  with  the  data  Lw  1  m.  With  the  position  of  this  knee  established,  7  w 
1.5  follows  from  matching  the  magnitude  of  cr^.  Variations  of  this  magnitude  in  the  wave- 
number  spectrum  over  short  ranges  are  consistent  with  the  observations  of  Briggs  [1989]. 

While  the  value  of  7  derived  from  Gensane’s  data  is  in  agreement  with  values  found  by 
Crowther  [1980],  the  characteristic  length  Lw  in  this  case  differs  greatly  from  the  typical 
value  of  100  m  found  by  Crowther.  This  difference  may  be  due,  in  part,  to  the  10  msec  pulse 
lengths  employed  in  Crowther’s  measurements  which  yield  insufficient  resolution  to  observe 
fluctuations  on  the  scale  of  a  meter.  On  the  other  hand,  in  Gensane’s  case,  the  very  limited 
range  scale  of  the  measurements  prevents  observation  of  large  patches. 

6.  Discussion 

A  theory  of  backscattering  from  slightly  rough  statistically  nonstationary  surfaces  has 
been  developed  and  applied  to  the  seafloor  reverberation  process.  The  validity  of  the  analy¬ 
sis  requires  only  that  the  dimensions  of  the  ensonified  area  be  large  on  an  acoustic  wavelength 
scale  and  that  the  surface  roughness  wave-number  spectrum  be  slowly  varying  near  the  Bragg 
wave  number  —  conditions  that  are  satisfied  for  seabed  reverberation.  Under  these  two  con¬ 
ditions,  the  scintillation  index  predicted  for  scattering  from  statistically  stationary  surfaces 
is  unity,  so  that  nonstationarity  must  be  invoked  to  obtain  agreement  with  experiment.  The 
theory  is  applicable  as  well  to  microwave  scattering  from  the  sea  surface  and  could  be  em¬ 
ployed  to  investigate  such  effects  as  the  influence  of  ripple  near  wave  crests  [Longuet-Higgins 
1987]  on  the  statistics  of  sea  clutter. 

When  the  theory  is  extended  to  the  time  domain,  another  feature  of  Bragg  scatter  be¬ 
comes  apparent:  that  the 'frequency  diversity  effect  observed  for  volume  scatter  and  scatter 


22 


from  very  rough  surfaces  at  high  grazing  angles  is  absent.  This  is  of  importance  for  seafloor 
reverberation  because  it  may  provide  a  method  of  distinguishing  the  contributions  of  rough 
surface  scatter  from  those  of  scattering  due  to  volume  inhomogeneities  within  the  seabed. 

In  this  analysis,  the  modified  reflection  coefficient  x(^)  of  (4)  has  been  treated  as  a 
constant.  However,  nonstationarity  of  this  parameter  could  contribute  significantly  to  the 
experimentally  observed  fluctuations.  A  theory  comparable  to  that  developed  here  that 
includes  variations  in  x  would  contain  a  strength  parameter  equal  to  the  normalized 
variance  of  |xp.  A  crude  estimate  of  this  variance  may  be  made  using  the  variances  of  sound 
speed  and  density  reported  by  Jackson  and  Briggs  [1992]  and  Yamamoto  [1995].  Assuming 
that  the  seabed  may  be  treated  as  a  homogeneous  halfspace  with  correlated  laterally  varying 
parameters,  for  sandy  sediments  is  estimated  to  be  on  the  order  of  unity  below  the  critical 
angle  and  significantly  lower  above  this  angle.  This  estimate  suggests  that,  if  sound  speed 
and  density  variability  is  a  contibuting  factor,  the  fluctuation  statistics  will  decrease  with 
increasing  grazing  angle.  Only  Crowther  [1980]  reports  such  a  dependence  of  fluctuation 
statistics  on  grazing  angle.  For  sands,  surficial  sediment  layering  [McDaniel  1992]  and  sound 
speed  gradients  [Moe  and  Jackson  1994]  may  strongly  effect  x  at  grazing  angles  above  the 
critical  angle.  Lateral  variation  of  these  parameters  would  hence  produce  an  increase  in 
fluctuation  with  grazing  angle.  It  is  evident  from  this  brief  discussion  that  careful  studies  are 
needed  to  fully  assess  the  contribution  of  sediment  variability  to  reverberation  fluctuations. 
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Figure  2.  Scintillation  index  for  scatter  from  statistically  stationary  sur¬ 
faces  having  Gaussian  and  power-law  wave-number  spectra. 
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Figure  3.  Ratio  of  the  scattered  intensity  to  that  predicted  at  the  Bragg  fre 
quency  for  scatter  from  a  stationary  rough  surface  characterized  by  a  power 
law  wave-number  spectrum. 
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Figure  4.  Scintillation  index  for  scattering  from  rough  surfaces  having  Gaus¬ 
sian  and  exponentially  correlated  nonstationarity  with  a  characteristic  length 
of  20  m. 


Figure  5.  Scintillation  index  for  scattering  from  rough  surfaces  having  Gaus¬ 
sian  and  exponentially  correlated  nonstationarity  with  a  characteristic  length 
of  5  m. 
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Figure  6.  Dependence  of  the  scintillation  index  on  beamwidth  for  scattering 
from  rough  surfaces  having  Gaussian  and  exponentially  correlated  nonsta- 
tionarity  with  characteristic  lengths  of  5  and  20  m. 
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Figure  7.  Normalized  temporal  covariance  of  the  impulse  response  for  scat¬ 
ter  from  a  rough  surface  having  power-law  wave-number  spectra  with  char¬ 
acteristic  lengths  of  20  and  100  m. 
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Figure  8.  Normalized  temporal  intensity  covariance  for  the  scatter  of  a  1 
msec  pure  tone  pulse  from  a  stationary  rough  surface  and  surfaces  having 
Gaussian  correlated  nonstationarity  with  characteristic  length  as  a  parame¬ 
ter. 
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Figure  9.  Normalized  temporal  intensity  covariance  for  the  scatter  of  a  10 
msec  pure  tone  pulse  from  a  stationary  rough  surface  and  surfaces  having 
Gaussian  correlated  nonstationarity  with  characteristic  length  as  a  parame¬ 
ter. 
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Figure  10,  Relative  predictions  of  <  p{ts  +  At/2)p*(rs  —  At/2)  >  by  (17) 
and  its  approximation  (18)  for  the  scatter  of  a  1  msec  pure  tone  pulse  from 
a  rough  surface  having  a  wave-number  spectrum  of  power-law  form  with  a 
characteristic  length  of  100  m. 
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Figure  11.  Comparison  of  measured  scintillation  indices  as  a  function  of 
effective  pulse  length  with  predictions  assuming  Gaussian  correlated  nonsta- 
tionarity. 
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